Anomalous quantum reflection as a quasi-dynamical damping effect 
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We develop a quasi-analytical theory for the quantum reflection amplitude of Bose-Einstein con- 
densates. We derive and calculate the decay-width of a Bose-Einstein condensate. A general relation 
between the time-dependent decay-law of the system and its quantum reflection amplitude allows 
us to explain the quantum reflection anomaly of Bose-Einstein condensates present in BEC-surface 
systems as a direct consequence of the repulsive particle interaction. 
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I. INTRODUCTION 



In recent experiments the quantum reflection of a 
sodium BEC from a silicon surface has been investigated, 
Surprisingly the quantum reflection probability fell 
off to zero at threshold and did not approach unity as ex- 
ected from corresponding single atom beam theory e.g. 
^ and experiments e.g., [1,H, 0, BQ- 



Theoretically so far, two numerical studies were able 
to qualitatively produce the quantum reflection anomaly. 
In[!l,i a model was used to motivate the quantum reflec- 
tion anomaly by introducing an empirical potential which 
modifies the atom-surface interaction. This lead qualita- 
tively to the observed effect. In the second numerical 
study an ab initio numerical simulation of the exper- 
imental situation reported in in three dimensions was 
used to explain the anomalous behaviour of the quantum 
reflection probability with the formation of scattering ha- 
los and the building of vortices that destroy the coherence 
of the condensate while being reflected. The anomalous 
behaviour of the reflection probability emerged when tak- 
ing into account only regions of the reflected BEC with 
more than 25% of its initial density. 



In the following we present a quasi-analytical approach 
to describe the anomalous behaviour of the quantum re- 
flection probability Due to its simplicity the ori- 
gin of the anomaly is elucidated. Furthermore, we iden- 
tify scattering along the axis of normal incidence on 
the surface as being mainly responsible for the anoma- 
lous behaviour of the quantum reflection probability. 
On base of this, we may use the framework of the 
well-understood spherical quantum reflection trap model 
pdl . [l^ . A spherical model has the advantage that there 
is only the normal direction for s- waves. This property 
strongly reduces the technical complexity of the general 
three-dimensional system and allows us to focus on the 
physics behind the quantum reflection of a BEC. Further- 
more, our model relies solely on the presence of particle- 
interaction and mean-field dynamics. 



II. QUANTUM REFLECTION 

For monochromatic atomic beams or linear wave- 
packets the universal law of quantum reflection is ob- 
tained by the fact that atom-surface potentials behave 
as step-potentials for incident momenta close to thresh- 
old, see d, 13, [ni, [l^. The universal threshold-law of 
quantum reflection is given by 



lim \R{k)\ = 1 - 2 6fc 

k 



(1) 



where the parameter b is called the threshold-length and 
k is the incident momentum. In case of an atom-surface 
potential it has been shown 0, 0, that close to thresh- 
old the retarded Casimir-Polder tail of the potential dom- 
inates the process of quantum reflection. The retarded 
Casimir-Polder tail is given by 
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where L denotes the location of the surface. For low 
incident momenta such potentials behave like 



U(z) = - — b-^t 



L), b = A 



(3) 



For the quantum reflection of linear wave-packets it has 
been shown in [ll| that the quantum reflection amplitude 
|i?(fc)| can be extended to the whole fc-space by 



\R{k)\ = 1 - 2 6fc « exp[-2bk] 



(4) 



with reasonable accuracy. It further has been shown that 
the decay-law for a linear wave-packet is linked to Eq. ([4]) 

by 



P{k,t) = \R{k)\^ 
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III. DERIVATION OF THE CONDENSATE'S 
DAMPING FUNCTION 

To study the influence of particle-interaction on the 
threshold-behaviour of quantum reflection, we turn to the 
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radial system that we have examined in detail from a dy- 
namical point of view, The radial Gross-Pitaevskii 
equation that describes the time evolution of an initial 
state ^(x, r = 0) in the atom-surface system in scaled 
form is given by 



■0(a;, t) 



' (x — 1) ip{x^ t) 



lV'(a;,T)| 
+ 7 5 ^{x,t) 



(6) 



Scaling is carried out by the intrinsic length of the system 
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In [12] we have shown that a repulsive self-interaction 
provides an additional contribution to the total energy 
of the system that is responsible for a faster decay of the 
particle-density than it is found in the linear case. Dy- 
namically, this is explained by a transformation of the 
self-interaction energy into kinetic energy as the initial 
state evoles in time. In fc-space, we interpret the pres- 
ence of the self-interaction as an additional barrier. The 
/c-modes close to threshold lie at the base of this bar- 
rier and are thus damped away by driving them over the 
edge of the potential. The effect of the self-interaction 
is strongest at the beginning of the time-evolution, be- 
cause its magnitude is directly related to the incident 
particle-density. For larger times, a considerable fraction 
of particle-density has already decayed, such that the in- 
fluence of the self-interaction becomes more and more 
negligible. 

For small times, we assume that the decay of the sys- 
tem can be described by a product-ansatz 

P(k,t;ct,7) = Pa{K,T) P^{k,t) = \R„{k) R^{k)\^'^'' 
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In Eq. 



the additional damping induced by the 



particle-interaction is described by Pj{fi,T). The 
product-ansatz of Eq. ([5]) reflects our assumption that 
for early times the effect of the step-potential and the 
particle-interaction can be treated independently. By 
this assumption, we derive the additional damping in- 
duced by the particle-interaction without taking into ac- 
count the presence of the atom-surface potential in Eq. 
®. 

The wave-function that solves Eq. ([5]) without po- 
tential can be decomposed into its momentum compo- 
nents by setting 

^(.,.)=/d.A(.)^.(.)expH^V]0,(.,r), 
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All effects of particle-interaction in Eq. ^ are described 
by (/). The Fourier-decomposition in Eq. ^ is defined by 



Jo 

and the basis of the system is 



Vk{x) ^ \I - smiKXl 



(10) 
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The decay can be calculated from the time-dependent 
dispersion, that is obtained by inserting the Fourier- 
decomposition Eq. ([9]) in Eq. ([6]), neglecting the atom- 
surface potential, multiplying by if from the left, in- 
tegrating out X and separating off the non-interacting 
parts. We thus obtain the time-dependent dispersion of 
the self-interacting system 



E^{k,t) = 7 / dKidK2 
Jo 



. dr (/>7(k, r) 
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(12) 



The function V{k, ki, K2) in Eq. (jl2p is a vertex-function, 
that describes the interaction of the modes of the system 
in momentum space. The main contribution to the vertex 
stems from the threshold region. We find 



V{k,Ki,K2) = 3k ^ A{ki) A{k2) - X 
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The decay of the system is described by the imaginary 
part of the dispersion, that can be identified by the line- 
width r.^. To calculate F-y, we expand E^{k^t)^ (/)^(k,t) 
in frequency-space, obtaining 

/d (jj 
- — Ej{k,uj) exp[— t] , 
2 TT 

/d uj 
— 0^(k,w) exp[-iw r] . (14) 



By using Eq. p4p we can perform a Laplace-transform 
on Eq. p2p . The imaginary part we are interested in can 
be extracted by using the well-known relation — = 
Pj^z^ — iTT5{Q, — Lo). We integrate out the (5-function 
and find 



7 / dni dn2V{n, Ki, K2) X 


— (/)*(AC2,W2)0-y(Kl,r2 - K2 +^2) (15) 

So far, we have only neglected the atom-surface potential. 
Since we are interested in the damping that is induced by 



the initial state we now make a relaxational ansatz that 
treats independent of lo by setting 



(16) 



We insert Eq. ([TB)) into Eq. IjlSp. carry out the remaining 
integration by the residual theorem and arrive at 
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Equation P7| is a self-consistent equation for the line- 
width of the self-interacting part of the system. From 
Eqs. p6l flT)) it follows, that the self-interacting part of 
the system initially decays like 



|0^(K,r)| K. 

Together with Eq, 
readily follows as 



P^(k,t) = exp[-2r^(K)T] . (18) 
. ([8]) the quantum reflection amplitude 



|i?(K;a,7)| = \R„{n)R^{n)\ 
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Following iJJj] the quantum reflection amplitude Eq. ^ 
describes the decay of the initial state in a uniform and 
quasi-stationary way by 



(p(r)} = / dn \A{n)\' |i?(«;; a, 7)!^ 
^0 



(20) 



where (p(r)) in Eq. ([20]) stands for the mean- value of the 
time-dependent particle density p(t) = dx |-!/'(a;, r)|^ 
that is obtained from numerics. 



IV. DISCUSSION 

Now we are going to test the fidelity of our theory by 
comparing it to numerical results. To make contact with 
our analysis in [Til . [T^ , we use the initial state 



vl'(a;, r = 0) = Mx exp [-ax]9[l - x] 



(21) 



where Af is the normalization constant and a is the dif- 
fuseness of the wave-packet. A diffuseness of a = 5, 
gives an initial kinetic energy i?kin = 1-5 x 10~^^ [a.u.] 
for sodium that corresponds to temperatures of approx- 
imatly InK, comparable to the experimental setup in 
[3,2']. As in 11] we chose a length-scale L = 4.47 x 
10^ [a.u.], and together with a potential-strength /34 — 
1.494 X 10" [a.u.] and mass m = 4.22 x 10^ [a.u.] for 
sodium we obtain a scaled potential-strength cr — 30. 
To illustrate the effect of particle-interaction, we chose 
7 = 0.5 . 
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FIG. 1: Damping-function r^(K) calculated by Eq. (I17I I for 7 
0.5. 



Figure ([T]) shows the additional damping that is in- 
duced by the particle interaction. It is clearly visible 
that the low-lying fc-modes are strongly damped, while 
higher /c-modes experience only a moderate to negligi- 
ble additional damping and thus show a regular decay 
according to universal quantum reflection Eqs. (UlIS]). 

Figure ^ shows a comparison between numerically 
calculated densities as function of the scaled time r and 
their mean-value approximation according to Eq. ([18]) 
for early times of the decay of the system. We have cho- 
sen a logarithmic ordinate for a better resolution. Fig- 
ure ^ shows clearly that after the einschwingvorgang 
is completed the mean-value approximation Eq. (|18[) in- 
terpolates the numerical curves. The approximation Eq. 
(|18p describes the decay in a quasi-stationary and uni- 
form way and thus does not reflect the oscillations of the 
numerically obtained data, that originate from the mo- 
tion of the wave-packet on the spatial range of the step. 
For times r > 0.35 our approach is not reliable anymore, 
because it takes only into account the damping accord- 
ing to the initial state, that naturally must overestimate 
the decay of the particle density for later times. A scaled 
time T ~ 0.35 as a range of fidelity corresponds to times 
t ~ 0.144 seconds. 

Figure ^ shows the quantum reflection probabilities 
calculated according to Eq. p^ . Figure ^ clearly 
demonstrates how the additional damping that is induced 
by the particle- interaction influences the quantum reflec- 
tion probability. The depletion of the quantum reflection 
probability for low-lying fc-modes corresponds to the re- 
sults reported in P, [3, [3 • The curves shown in Fig. ([3]) 
must be interpreted as a quasi-stationary mean-value of 
the quantum reflection probability that is to be expected 
for times r < 0.35. For later times, as the influence of 
particle-interaction declines along with the decaying den- 
sity, we expect a reduction of the anomalous behaviour 
that, for very large times, will experience a cross-over to 
the universal law of quantum reflection Eqs. (HKH]). 
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FIG. 2: Comparison between numerically obtained decaying den- 
sities p(t) (full lines) and approximated decaying densities (p{t)) 
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FIG. 3: Quantum reflection probabilities according to Eq. (|19|l for 
(cr = 30, 7 = 0.5) (full line) and (cr = 30, 7 = 0) (dashed line). 



V. SUMMARY AND CONCLUSION 

Together with our investigations on the dynamical 
properties of a BEC-surface system, fT5], the behaviour 



of Eq. (|19p as shown in Fig. Q proves, that quantum 
reflection has only a limited effect on BECs with repul- 
sive particle-interaction. We thus explain the quantum 
reflection anomaly by the fact, that especially the low en- 
ergetic, collective components of the wave-packet moving 
normal to the surface react extremely sensitively on the 
presence of an additional positive energetic contribution 
as a repulsive particle-interaction provides. The crucial 
role of the repulsive particle-interaction and its transfor- 
mation into kinetic energy was already mentioned in [l[ 
and theoretically demonstrated in The influence of 
the repulsive interaction on higher energetic components 
is by far less dramatic, explaining the regular behaviour 
of Eq. (fT9|) for higher momenta. The fact that we can ex- 
plain the quantum reflection anomaly with the behaviour 
of the line- width as function of momentum may give hints 
for future experiments to measure the line-width of the 
condensate instead of the reflection-probability, which, as 
emphazised in [1, 2], is quite sensitive to external influ- 
ences. 

The purpose of our present theory is neither to exactly 
recalculate the experimental data from [l], nor the 
numerical results reported in [l^, but to elucidate the 
physical mechanism that is responsible for the anoma- 
lous behaviour of the quantum reflection probability on 
the basis of a simple but well-understood model. Our 
quasi-analytical theory focuses exclusively on the direc- 
tion of normal incidence, which we have assumed and 
proved to contribute signiflcantly to the quantum re- 
flection anomaly. Our analysis demonstrates that the 
anomaly exists and that it may be found in experiments. 
Furthermore, the experiments reported in [J Q so far 
contain enough uncertainties to doubt about the rele- 
vance of a direct and quantitative comparison between 
theory and experiment in the present state of the art. 
However, our analysis may encourage further experimen- 
tal work in this fleld. 

In the framework of the simple spherical model, p^ . 
we have developed a quasi-analytical theory for the 
anomalous behaviour of the quantum reflection proba- 
bility. Our theory explains the anomalous behaviour of 
the quantum reflection probability of BECs as a direct 
consequence of the repulsive particle-interaction. The 
key-quantity for understanding the quantum reflection 
anomaly of BECs was shown to be the line-width of the 
condensate as function of momentum. 
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